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Abstract 


An improved method for digital interpolation of signals for a second interpolation filter is disclosed which 
permits a high signal/noise ratio with a minimum amount of circuitry for an overall system comprising first 
and second interpolation filters. The method for digital interpolation of signals requires multiplying delayed 
input values locked to a first signal by corresponding weighting factors which are dependent on a time- 
difference value determined by the interpolating instant and the time grid of the first clock signal. The 
weighting factors are determined by an impulse response in the time domain. The associated transfer 
function has an attenuation characteristic in the frequency domain which, with respect to the stop bands, is 
limited essentially to the alias regions located at the frequency multiples of the first clock signal. Each of 
these alias regions is assigned at least two adjacent zeros, or in the presence of double-order zeros, at least 
one of the alias regions and the associated periodic alias regions are assigned at least one further zero of 
the transfer function. 
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(57) Verfahren zur digitalen Interpolation von Signa- 
len durch Multipiikatioh verzOgerter, mit einem ersten 
Taktsignal (f1) verkoppelter Eingangswerte (x(k)) mit 
entsprechenden Gewtchtungsfaktoren (Ci, C|). die von 
einem Zeitdifferenzwert (At; dt) abhdngig sind, derdurch 
den Interpolationszeitpunkt (tP) und das Zeitraster des 
ersten Taktsignals (f 1) bestimmt wird. wobei die Gewich- 
tungslaktoreh (01; G|) durch eine Impulsantwort h(t) im 
Zeitbereich bestimmt werden. Die zugelidrige Obertra- 
gungsf unktion H(f) weist im Frequenzbereicli ein Signal- 
ddmpfungsverhalten auf, ' das bezuglich der 



Sperrbereiche im wesentlichen auf die bei den Fre- 
quenzvielfachen des ersten Taktsignals (ft) liegenden 
StOrbereiche (1S. 2S, ...pS. ...) beschrankl ist. Jedem 
dieser Stdrbereiclie werden mindestens zwei nebenein- 
anderliegende Nullstellen zugeordnet, Oder beim Vpr- 
handensein von Nullstellen doppelter Ordnuhg wird 
mindestens einem^der Stdrbereiche (pS) und den zuge- 
fidrigen periodischen Stfirbereichen (3pS. 5pS, 7pS, 
9pS. ...) mindenstens eine wettere Nullstelle der Ober- 
tragungsfunktion H(f) zugeordnet. 
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Beschreibung 

Die Erfindung betrifft ein Vafahren zur digitalen Interpolation von Signalen durch Multiplikation verzogerter, mit 
. einem ersteh Taktsignal verkcppelter digitaler Erngangswerte mit entsprechenden Gewichtungsfaktoren (= Filterkoeffi^ 
5 z ienten) und Aufaddierung der multiplizierten Werte, wobei die Gewichtungsfaktoren von einem Zeitdiff erenzwert abhSn- 
. gig sind, der durch den Interpblationszeitpuhkt und das Zeitraster des ersten Taktsignals bestimnit wird. vergleiche den 
Oberbegriff des Anspruchs 1 . 

Derartige Verfahren werden bei Interpolations! iltern venA^endet. die insbesondere als Teilschaltungen in digitaJen 
Schaltungssystemen dienen, bei denen eine betiebige Anderung der Abtastrate von digitafen Signalen erfbrderlich ist. 
10 Systeme, die sich hur mit einfachen ganzzahiigen Abtastratenverhaltnissen befassen, sirvJ nicht Gegenstand der Erfin- 
dung. Bei beliebigen: Abtastratenverhaltnissen soil das mit einem ersten Takt abgetastete Signal in eine Abtastfblge 
umgewandeltwerdea die mit einem beliebigen zweiten Takt verknupfl ist, der im ungunstigsten Fall sogar ein irratiqnales 
Frequenzverhdltnis zum ersten Takt aufweist. Auch sollen mit der Erfindung diejenigen Zahlenverhaltnisse realisiertsar 
sein. bei denen das Meinste gemelnsame Frequenzvieifache beider Takle bei einer unenvdnscht hoheh Frequenz liegt. 
15 Aus "IEEE. Transactions of Acoustics. Speech and Signal Processing", Band ASSP-32. Nr. 3, Juli 1 984, Seiten 577 
bis 591 sind unter dem Titel "Digital Methods for (inversion Between Arbitrary Sampling Frequendes", Alitor: T A. 
Ramstad, Verfahren zur beliebigen Anderung von Abtastraten beschrieben. Die zugehOrigen Sishaltungen weiden als 
hybridiB Systeme bezeichnet. die aus einem ersten Interpolationsf ilter mit festem Abtastratenverhaftnis und einem zwei- 
ten Interpolationsf liter (= ahaloger Resampler) bestehen. Mitdem zweiten Interpolationsf ilter werden die Zwischenwerte 
20 bestjmmt, die zeitiich beliebig zwischen den festen Abtastwerten des Abtastrasters hach dem zweiten Interpplationsfilter 
liegen uiid damit beliebige Abtastfatenyerhgltnisse zulassen. Das erste Interpolationsf ilter enthalt als.Kombina^^^ eine 
. lirrterpolatiohseinrichtuhg und ein djgitales Filter; Mit der Intei^piDlatioris^iiiricht^ die auch ate UbeVabteisteinrichtung 
bezbichnetwird';.werden "Nuir-W^ einem Uberabtastfaktqr N zwischen die ursprdnglichenAb^ 

eingefugt; Mit eiher Dezimiereinrichtung werden gegebenenfalls ursprungliche Abtastwerte unterdruckt. Erst das nach- 
25 folgende digitale Filter giattet in beiden Fallen den Verjauf der digitalen Abtastwerte. wobei insbesondare de SprQnge 
auf die Null-yVerte ausgegOchen werden, so da6 das Spektrum des fslutzsignals njcht durch hdhere frequenzkompo- 
nenten verfSlscht wird. Hierzu ist das erste Interpolationsfilter so ausgebildet. da3 grdlBere Frequenzbereichsluckeri tn 
dem sich bis ins Unendliche erstreckende Spektrum gebildet werden. Denn auch bei der Oberabtastung gilt, daB sich 
die Freqiienzspektren bei der halben ursprunglichen Abtastfrequenzen und.deren Vielfachen spiegeln. Nach der Inter- 
30 poiationseinrichtung und pach dem digitajen Filter ist allerdings von einer heuen Abtastfrequenz - im Folgenden auch 
erster Abtasttakt genannt - auszugehen. die irn ganzzahiigen Frequenzverhaithis zui' usprunglichen Abtastfrequenz 
. steht. Das digitale RIter entfernt dabei wie ang^geben die verbleibenden spektralen Kbmponenten zwischen dem Nuti- 
signalband und dem gespiegelten Band bei der neuen Abtastfrequenz und den zugehOrigen FrequenzvieHachwerten. 
Die Funidionsweise des digitalen RHers ist dabei eihfach diejenige eines digitalen Tiefpasses, der den Nutzsignalbereich 
35 durchldBt und die daruberliegenden Frequenzkomponenten unterdruckt, wobei allerdings entsprechend d^ni Abtast- 
theorem die bekannte Spiegelung bei der halben Abtastfrequenz auftritt. Ein digitaler Tiefpa3 l^n also niemals die 
Vielfachen der Abastfreqenz unterdrucken. 

Die spektralen Kbmpphenten b^i der neuen Abtastfrequenz und den, Frequenzvielfachwerten mussen fur die Rea- 
lisierung.beliebiger Abtastratenyerfiaitnisse durch das zweite. an^loge Interpolationsfilter unterdruckt yverden. Werden 
40 diese Stdrkqmpohehteh hicht uhterdrQckt/dann werden bieider Erzeugung beliebiger Abtastratenve^^ 

ponenten im Nutzsignalband erzeugt. Die Theorle des ersteh Interpolatiqrisf liters ist hihreichend in "Proceedings of The 
IEEE**, Band 61. Nr. 6, Juni 1973, Seiten 692 bis 702 in dern Aufsatz'"A Digital Signal Preceding Approach to Interpo- 
lation" von R W. Schafer und L. R. Rabiner beschrieben. 

FQr das zweite Interpolationsfilter gibt T. A. Ramstad in dem oben genannten Aufsatz in "IEEE, Transactions Oh 
45 ' Acoustics. Speech and Signal Processing" . Juni 1984, Seiten 577 bis59l folgende MSglichkeiten an: 

'1 . Eine Abtastiialteschaltung - dieses Filter hat einfache Nullstellen in der Mitte der noch verbleibenden spektralen 
Komponenten bei den der heuen Abtastfrequenz. Diisser Filtertyp wird auch ats einfaches Kammfilter bezeichnet. 
2, Eine lineare Interpolationsschaltung - dieses Filter hat doppelte Nullstellen in der Mitte der noch verbleibenden 

50 spektrialen Komponenten bei den Frequenzvielfachwerten der neuen Abtastfrequenz. Damit wird eine wirksamere 
UnterdnCtckung der Steranteile als bei der einfachen Abtasthalteschaltung erreicht. Diese erfordert beispielsweise 
fur ein Signai/Rauschverhg^ltnis von 1 0OdB einen Uberabtastfaktor N von 200, der durch das erste Interpolationsfilter . 
realisiert werden muB. In der genannten Uteraturstelle wird ein "optimales" Interpolationsfiltei^ angegeben, das eine 
Absenkung des Oberabtastfaktors auf einen Wert von 164 ermoglicht. Selbst die Venwendung sehr rechenaufwen- 

55 digerrLagrange-lnterpolationsfilter hoherer Ordnung eriaubt keine wesentiiche Absenkung der erforderlichen Ober- 
abtastung. . Sowohl das lineare Interpolationsfilter als auch das Lagrangerlnterpolationsfilter zeigen ein 
Karhmfilterverhalteh. 



2 



EP 0 696 848 A1 



Aus EP-A 0 561 067 1st ebenfalls ein Verfahren mit einem hybriden System zur Abtastratenumsetzung bekannt. 
Dieses System arbeitet allerdings nur mit einem Uberabtastfaktor N = 2 und erreicht damit nur einen relativ schlechtes 
Stgnal/Rauschverhaitnls. Dies ist jedoch zulSssig, da dieses System fOr Videosignalanwendungen gedacht ist. Das 
zweite Interpolationsfilter ist dabei als Tiefpa3 realisiert. der alle Frequenzen oberhalb dem 1,5-fachen Wert der 

5 ursprunglichen Abtastfrequenz unterdruckt. Das analoge TiefpaRverhalten wird mit einem Transversalfilter erreicht, bei 
dem die Gewichtungsfaktoren der gespeicherten Abtastwerte von einem Zeitdiffererzwert abhSngig sind. der im folgen- 
. den auch als Inter-Sample-Position bezeichnet wird. Ein derartiger TiefpaB unterdruckt dabei nicht nui^ die verbleibenden 
spektralen Komponenten bei den Frequenzvielfachwerten der neuen Abtastfrequenz sondern den gesamten Spektral- 
bereich oberhalb seiner Sperrflanke. Bei einem vergleichbaren DurchlaB/Sperrverhalten ist ein derartiger TiefpaB immer 

10 aufwendiger zu realisieren als eine entsprechende Kammfilteranordnung. Der TiefpaB yerarbeitet zwardigitale, abge- 
tastete Werte - In seiner Wirkung stelH er jedoch einen analogen TiefpaB. einen analogen Resampler, dan 

Aus "Journal of Audio Engineering Society*. Band 41, Nr. 7/8. 1993. Seiten 539 bis 555 ist von den Autoren R. 
Adams und T. Kwan unter dem Titel "Theory and VLSI Architectures for Asynchronous Sample-Rate Converters" ein 
Verfahren fOr ein Abtastratenwandlungssystem bel^nnt, das einerseits die Ven/vendung einfacherer Abtasthalteschal- 

IS tungen und andererseits die Verwendung von TiefpaBfiltern als analoge Resampler behandelt. Fur einfache Systeme 
werden wegen des geringen Schaltungsaufwandes Abtasthalteschaltungen empfohlen - f Qr hdherwertige Systeme wird 
dagegen die Ven/vendung eihes TlefpiaBf liters nach dem ersten Interpdiationsfilteir vorgeschlagen, das den gesamten 
Frequenzbereich oberhalb des Nutzsignalbereichs unterdrOckt. 

Aus den bekannten Literaturstellen-geht hervor, daB nach der N-fachen Oberabtastung urid Filterung nach diem . 
' 20 ersten Interpolationsfflter im Frequenzspektrum auf jeden Fall StSrbereiche vorhanden sind, deren Mittenfrequenzen 
bei .deri Frequenzvielfachwerten der neuen Abtastfrequenz liegen,; Die Bandbreite jedes StOrt>ereichs ist dabei gieich 
derdbppyten' BiandbriBite des Nutzsignals/ Wehn die Nyqui^tbed pigitalisieruhg erfQIlt ist, 

dann hat die Biandbreite des.StOrbereichs im Grenirfkll'msodmal den W^rt der ursprunglichen 'A^^ tige 
und Bandbreite sdmtliche'r StOibereiche ist im Frequenzspektrum durch die ursprungiiche Abtastfrequenz und den Ober- 

25 abtastfaktor N definiert. Die N-fache Oberabtastung der ursprunglichen Abtastfolge bewirkt. daB die relative Breite der 
Stdrbereichebezogenauf die neue Abtastfrequenz umden Faktor 1/Nreduziertwird. Di€|serleichtertdieTrennung des 
Nutzsignalbandes vom jeweiligen StOrbereich, well der Obergangsbereich zwiischen DurchlaB- und Spenrbereich fur das 
zweite Interpolationsfilter vergrdBert wird, wodurch der erforderliche Schaltungsaufwahd fur das zweite Interpolations- 
filter verringert wird. Dies wird jedoch durch einen hOheren Aufwand fOr das GIdttungsf liter im ersten Interpolationsfilter 

30 erkauft. 

In der Regel hat man entweder ein sehr aufwendiges erstes Interpolationsfilter, z. B. durch eine hohe Oberabtastung. 
und ein einfaches zweite Interpoiationsfilter, z.B. einen linearen Interpolator, oder aber man hat ein einfaches. erstes 
Interpolationsfilter. z.B; mit einer geringen oder gar keinen Oberabtastung. und ein sehr.aufwendiges TiefjDaBfilter mit 
dem der analoge Resampler realisiert wird.; In der letzgenannten Literaturstelle werden beispieisweise 16.000 feste 

.35 Stutzstellen fur das aufwendige TiefpaBf liter gespeichert 

Aufgabe der Erf indung ist es daher, fur das zweite Interpolationsfilter ein verbessertes Verfahren zur digitalen Inter- 
polation von Signalen anzugeben, welches eih hohes Signal/RauschverhdItnis ermOglicht. Dies soli mit einem mOgfichst 
geringen Schaltungsaufwand fur das Gesarntsystem aus erstem und zweitem Interpolationsfilter realisiert werden, 
Die Ldsung der Aufgabe ergibt sich aus den Merkmalen des Anspruchs 1 . Vorteilhafte Weiterbildungen der Erfin- 

40 dungen ergeben sich aus den UnteransprOchen. 

Die Erf indung bestehtdarin. daB sich die Signalddmpfung des zweiten Interpolationsfilters im wesentlichen kamm- 
filterartig nur auf die Stdrbereiche bei den Fr^uenzvielfachwerten der neuen Abtastfrequenz beschrdnkt und fur die 
dazwischenliegenden Bereiche einen beliebigen Signalverlauf zul3Bt. Dabei weist die Obertriagungsfunktion des zweiten 
Interpolationsfilters bezuglich jedem der Stdrbereiche miridestens zwei Nullstellen auf, wobei die mindestens zwei Niillsi- 

45 ellen in jedem Stdrbereich nebeneinander liegen oder falls es sich bei den mindestens zwei Nullstellen um elrie NuUstelle 
zweiter Ordnung (= zweiten Grades) handelt, mindestens einem der Stdrbereiche mindestens eine weitere Nullstelle 
der Obertragungsfunktion zugeordnet wird. Diese weitere Nullstelle wird vorzugsweise, zurgezielten Verbreiterung des 
Sperrbereiches. neben die vorhandene Nullstelle zweiter Ordnung gesetzt. Durch die mindestens eine weitere Nullstelle 
werden ferner bei den zugehdrigen Frequenzvielfachwerten weitere Nullstellen erzeugt, die dann in der Regel innerhalb 

50 der zugehGrigen periodlschen Stdrbereiche liegen. Dadurch wird die UnterdrQckung der unenwunschten Frequenzbe- 
reich weiter verbessert. Zudem nimmt mit steigenden Frequenzen auch in den DurchlaBbereichen des Kammf liters 
dessen Ddmpfungsverhalten zu. so daB die Lage der Nullstellen innerhalb der Stdrbereiche dort immer unkritischer 
wird. Durch die jeweilige Obertragungsfunktion werden die Gewichtungsfaktoren der miteinander zu verknOpfenden 
Abtastwerte im zweiten Interpolationsfilter bestimmt. wobei die Gewichtungssfaktoren wiederum von dem jeweiligen 

55 Zeitdifferenzwert des Abtastzeitpunktes zwischen dem ersten und zweiten Abtastraster abhangig sind. 

- Die welteren Nullstellen oder die Modif ikation Ihrer Position lassen sich durch folgende Verfahren erreichen, wobei 
die Verfahren miteinander kombinierbcU^ und aitich mehrfach aufeinanderfolgend anwendbar sind: 

a) die Verwendung einer Abtasthalteschaltung hdherer Ordnung - also mit mindestens dreifachen Nullstellen, 



3 



EP 0 696 848 A1 



b) die Vefschiebung und Addition der Impulsantwort im Zeitbereich und Bildung einer resultierenden, neuen Impuls- 
antwort und schlleBlich 

c) die Verschiebung und Muttiplikation der Ubertragungsfunktion im Frequenzbereich und Bildung einer resultieren- 
den, neuen Impulsantwort, wodurch in jedem Sperrbereich mindesten zwei nebeneinanderliegende Nullstellen 
erzeugt werdeh. 

Die Verfahren nach a) und b) lassen sich mittels einer besonders geeigneten Filterstrukiur, der als Name die Bezel* 
chiing "Switchable Time-Continuous FIR" (= STC-FIR) gegeben wird, realisieren. 

Der Zeitdiff erehzwert bzw. die Ihter-Sample-Position definiert einen beliebigen Zeilpunkt zwischen zwei festen Abt- 
astzeitpunkten des ersten Frequenzrasters: Durch das erflndungsgemSGe Verfahren laBt sich die weitere bzw. lassen 
sich die werteren Nulisteilen bezuglich der Storbereiche beliebig positionieren, so daB eine gezielte St6runterdruckung 
mdglichwird. Durch einienseitlichen Versatz der Nullstellen imjeweiligen StdrbereichlaBtsich der jeweilige DSmpfungs- 
bereich einfach in der Breite modifizieren. Durch gezielte Positionierung^weiterer Nullstellen kann man beispielswei^e 
den Obei^abtastungsfaktor im ersten Interpolatibnsf liter von N = 200 auf N s 8 reduziereii. wobei das zwerte Interpola- 
tionsf ilter nur 5 Abtastwerte miteihahder verknupfen muB. Unter cGesen Bedingungen laBt sich ein Signal/Starverhaitnis 
von lOOdB erreichen. Dadurch wird der notwendige Schaltungs- bzw. Programmieraufwand fbr das Abtastratehr 
Umwandlungssystern stark reduziert und die Taktfrequenzen bleiben zudem niedrig. 

Die Erfihduhg und.vbrteilhafte Ausgestaltungen werdeh nun anhand der Figuren'der Zeichhung ndtier erldutert': 

Fig. .1 zeigt schertiatisqh die grundsatzjiche Funktiohsweise von Verfahren zur beliebigen Abtastratenurhseteijng. 

'.Fig^S v^ransc^^ J V ' ' . ' 

Fig. 3 zeigtjschenr^ . :. ; » 

Fig. 4 ?eigt die IrnpulsantWdrten.j^^^^ ersteir, zweiter, dritteir und:yiei1eti^ oWijrig^'.gm dem- 

ersten Verfahren a)/ , • ' " : * v 

'Fig. 5 zeigt die entsprechenden Ubertragungsfunklionen zu den Abtastfunktibnen hSherer Ordniingi 
Fig; 6 zeigt die Impujsantworten gem^B dem zweiten. Verfahren b). das eine Verschiebung im Zeitbereich betrifft, 
Fig. 7 zeigt die entsfjredienden Obertraguhgsfunktionen zu den Zeitverschiebungen ta = T/2 . T/4'und T/8, 
Fig. 8 zeigt ein Beispiel fur eine resultierende Obertragungsfuhktion gemaB dem dritten Verfahren c), das- eine- 
Verschiebung im Frequenzbereich betrifft, 

Fig. 9 ziBigt die entsprech'ende lrripulsant>yort zu der Frequenz\^^ . 

Fig . 1 0 zeigt schematisch als BlocksqhaHbild eiri^e Realisierung des zweiten Interpolatibnsfilters dhnlich einem FIR- 
" Filter. ^ ' ^ ^ ■ . . V'- 

Fig.. 11 zeigt schematisch als Blockschaltbild ein vorteilhaftes AusfQhi^ijngst)eispiel des tnterpolationsf liters mit einer 
umschaitbaren Matrix, * .. 

Fig. 12 zeigt schematisch im>Zeitdiagramm diie OberfOhriing von einem ersten Abtastraster in ein beliebtges zweites 
Abtastraster und * ; - 

Fig. 13 zeigt schematisch wie in einer Matrix die verschiedenen Gewichtungsfaktoren in Abhdngigkeit vpn einem 
Ze^tdifferenzwert gebildetwerden. ; . . • f 

In Fig: 1 ist als Blocksctialtbild ein. Scherna fdr beliebig Abtastratehumsetzungen dargesteltt/wie es Qbiicherweise 
zur Eriauterung des^Veffahrens verwendet wird: Ein ahaloges Eingarigssignal Sej^ wird mittels 'eines A^^ 
wandlers 1 digitalisiert, wobei die Digitalisierungsfrequehz f^in istl Seibstverstdndiich muB die Nyquist^ erfulit 
sein. so daB sich der Frequenzbereich des analpgen Signals Sejn hdcHstens'bis zur halben Abtiastfrequenz fej^/S 
erstreckt. Nach diem Analog-Digitalwandler tomnit eine Interpolatiohsstufe 2, welche die ursprQngliche Abtastrate fgi^ 
um den Faktor N erhOht. Hierzu werdeh in der Regel "Nuir-Werte zu den zus§tzlichen Abtastzeitpunkten als Abtastwerte 
eingefugt/ Durch die Digitalisierung mit der Frequenz fein spiegelt sich das Nutzsignalspektrum bei der halben Abtast- 
frequenz fejn/2 urid wird zwischen deri Frequenzwerten f^j^ und 1,5'fein wiederholt und so fort, vgl. Fig. 3. Durch die 
Uberabtastung in einemjnterpolator 2.1 mit dem ganzzahligen' Faktor N §ndert sich an dieser pertodischeii Folge des 
Signalspektrums nichts. Ein digitales Rlter R entfernt jedoch aus den Storspektren moglichst alie Frequenzarrteile, die 
zwischen der halben ursprCinglichen Abtastfrequenz und dem ersten Storbereich bei der neuen Abtastfrequenz N'fej^. 
und zwischen den weiteren Stbrbereichen liegen. Damit liegt entsprechend Fig. 3 ein Signal vor, das nur den Frequenz- 
bereich des Nutzsignals enthdit und daruber hinaus nur Storahteile, die bei den Vielfachwerten der neuen Abta^requenz 
liegen. Wenn dieses Signal mit einem Digital-Analogwandler 3, der mit der Abtastfrequenz f 1=N-f , arbeitet. zurQck- 
gewandelt wird, dann liegt nach einer entsprechenden Giattung angenShert der ursprungliche SIgnalverlauf Sejn'wieder 
vor. Der Digital/Analog- Wandler 3 erzeugt dabei theoretisch nur zu den unendlich kurzen Abtastzeitpunkten der Frequenz 
f1 Analogwerte x(t), die durch ein Giattungsfilter h in ein kontinulerliches Signal y(t) umgewandelt werdeh. Das Glat- 
tungsf ilter h. das wesentlicher Teil des zweiten tnterpolationsf liters. 4 ist, kann daher zutreffend audi analoges tnterpo- 
lationsf ilter genannt werden. 
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Fur das Verstandnis der Abtastratenumwandlung ist davon auszugehen, da3 das Signal y(t) nach dem Glattungs- 
filter h mdglichst dem Eingangssignal Sein gleich ist. das heiBt da6 alle Nutzsignalkomponenten vorhanden sind und 
keine StOrsignaikomponehten. Ein unterschiedlicher Skalierungsfaktor ist dabei unerheblich. Wenn diese Bedingungen 
ideal erfulltwaren, danri kOnnte mittels einer zweiten Analog-Digital-Wandlereinrichtung 5 das kontinuierliche Signal y(t) 

5 mit einem beliebigen zweiten Abtasttakt f2 digitalisiert werden. Die digitalisierten Werte wurden eine zweite digitalisierte 
AjDtastfolge y(d) lief em. Wenn das Signal y(d) spdter hriit einer niederen Ablastfrequenz fgus weiter verarbeitet werden 
soil, dann schlieBt sich der Digitalisierungsstufe ein digitales TiefpaBfilter Fo und eine Dezimierstufe 6 an, wodurch der 
Abtasttaktf2 umdasganzzahlige Freqiienzverhdltnis No reduziert wird. Dem Ausgangssignal Saus ist die Abtastfrequenz 
faus zugeordnet, die sich aus der zweiten Abtastfrequenz f2 mit = f2/No , ergibt. 

10 Aus Fig. 1 ist ohne weiteres erkennbar, daB die Genauigkeit der Abtastratenumwandlung im wesentlichen von der 
Gehauigkett des Signals y(t) abhdngig 1st, und dieses Signal hSngt wiederum von den Eigenschaften des Giattungsf ilters 
h und des digitaleh Filters Fi ab. Bet der praktischen Realisierung verzichtet man selbstverstdndlich auf den analogen 
Resampler R in der dargestellten Form, weil er eine zusdtzliche Digital/Analog- und Analog/Digital-Umsetzung fur das 
analogie GIdttungsfilter h erfbrdert. Es sind ndrnlich, wie eingangs beschrieben, auch analoge GIdttungsverlahren 

IS bekannt, die mit einzelnen Abtastwerten auskommen und die erforderlichen Zwischenwerte durch Interpolationverlahren 
erzeugen. Den einfachsten Fail stellt hierbei die einfache Abtasthalteschaltung (= latch) dar. die den letzten Abtastwert 
der Folge x(k) so lange festhSlt. bis der neue Abtastwert vorliegt. Eine A/erbesserung bringt der llneare Interpolator, der 
die beiden letzten Abtastwerte speichert und daraiis durch liheare Interpolation einen Zwischenwert zum gewOnschten, 
dazwischenliegenden Ihterpolationszeitpunkt tP blldet. Andere bekannte Gl§ttungsfilter sind sehr auwendige TiefpaB- 

20 filter Oder das Lagrange-Filter Da bei alien diesen Filtern die Inter-Sample-Position beliebrg w^hlbar ist wird digital auf 
diese Weise ein analoges Rlter realisiert, daher auch die Bezeichnung "analoger Resanripler". \^ . - 

In Fig. 2 ist yereinfacht das kontinuierliche Signal y(t) uber der Zeitachse t zusanrin^ien mit dem,erstai und zweiten 
Abtastraster f i bzW. f2 dargestelit. Es soil beispielsweise fur deh Zieitpunkt tP eine neuer Abtastwert P2 gebildet werden. 
Mit dem analogen Resampler R nach Fig. 1 erfolgt dies einfaich durch den zweiten Analog/Digital- Wahdier 5, der vom 

25 zweiten Abtasttakt f2 zum Zeitpunkt tP getaklel ist. Ohne den Umweg uber den analogen Resarhpler R muB der Abt- 
astwert P2 jedoch aus benachbarten Abtastwerten P11. PI 2, P13, P14 der ersten Abtastfolge interpoliert und damit 
berechnetwerden, Derzeitlich unmittelbarvorausgehende Abtastwert P12 uhdder InterpolationszeitpunkttP bestirhmeh 
dabei einen Zeitdiffeirenzwert At (in der folgenden Beschreibung wird hierfur meist die Bezeichnung "dt" verwendet ) , 
der zur Berechnuhg des interpotierten Abtastwertes P2 dient. Bei der linearen Interpolation bestimmt er ndmlich den 

30 Anteilsfaktor aus den beiden benachbarten Abtastwerten P12 und P13. Fur die lihearen Interpolation sind die weiteren 
Abtastwerte P1 1 und PI 4 und alte weiteren Werte ohne Interesse. Der Zeitdifferenzwert dt bezieht sich selbstverstdnd- 
lich auf die PerioderKlauer T1 des ersten Abtasttaktes f1. In der Zeitdarsteliung stellt sich die lineare Interpolation als 
symmetrische, dreiecksfdrmige Impulsantwort h2(t) im Zeitbereioh dar, wobei die Symmetrieachse durch den Interpo- 
lationszeitpunkt tP definjert ist. Die zeitliche Breite der impulsantwort h2(t) umfaBt dabei zwei Abtastperioden T1 und 
damit:zwei gespeicherte Abtastwerte. Der Anteil- Oder Gewichtungsfaktor ergibt sich aiis dem Schnittpunkt der Abtast- 
werte P12 und P13 mit der Impulsantwort h2(t}. In Fig. 2 fdllt der Abtastpunkt P1 4 es ersten Abtasttaktes f1 mit einem 
Abtastpunkt des zweiten Abtasttaktes f 2 zusammen, so daB dort der Zeitdifferenzwert dt zii Null wird. Die Berechnung 
des jeweiiigen Ahteilfaktors ist In alien Fallen unabhangig von der Fr^quenz des zweiten Abtasttaktes f2. Lediglich der 
Interpolationszeitpunkt tP isl fOr die Berechnung maBgebend. 

40 \n Fig, 3 sind schematisch eiriige zu Rg. 1 gehfirende Frequenzspektren dargestelit Das Nutzsignal Kat entspre- 
chende dem Nyquist-Kriterium hochstens die Bandbreite B=f ^j^/ 2 . Dem entspricht eine doppelt so groBe Bandbreite 

Bs bei den Storbereichen I S pS, die mittig zu den Frequenzvielfachwerteh des ersten Abtasttaktes f1 liegen. In der 

ersten Zeile sind gestrichelt die Frequenzspektren angedeutet, die zwischen dem Nutzsignalband und den jeweiiigen 
Stdrbereichen liegen und die infolge der Spiegelung an der halben ursprunglichen Abtastfrequenz fejn/2 entstehen. 

45 Durch die Oberat>tastung mit dem Faktor N dndert sich dabei nichts. Die gestrlchelteh Frequenzberelche sollen jedoch 
mdglichst gut durch das digitale Filter Fi unterdruckt werden. Das digitale Filter Fi ist damit ein typlscher digitaler TIelpaB, 
dessen Frequenz sich bei der halben Abtastfrequenz H/Z spiegelt. 

In der zweiten Zeile stellt Fig. 3 den Frequenzverlauf und damit die Ubertragungsfunktion H(f) im Frequenzberelch 
eines mdgtichst idealen • und darhit analogen - Tielpasses dar, der als Interpolationsfilter 4 alle Stdrsighale oberhalb 

so des Nutzsignalband es B unterdruckt. In der darunterllegenden Zeile wird der gewunschte Frequenzverlauf H'(f) eines 
GISttungsf liters h dargestelit, bei dem kammf ilterartig nach der Erfindung im wesentlichen nur die periodischen Stdrbe- 
reiche IS, 2S, pS unterdruckt werden. Der D^mpfungsverlauf in den dazwischenliegenden Berelchen ist uner- 
heblich, wenn sie vom digitalen Filter Fi vorher genugend unterdrQckt wordeh sind. In der letzten Zeile wird schennatisch 
die spektrale Verteilung y(f) des interpolierten Signals y(t) dargestelit, die bei den einzelnen Stdrbereichen nur noch 

55 kteine Signalreste enthdit, die kleiner als das gewunschte Signal/Rauschverhaltnis sein sollen. Anstatt der Bezeichung 
H(f) wird insbesondere in den nachfolgenden Gleichungen die gleichwertige Bezeichnung H(g>) verwendet mit 

<0 = 2- 7C-f. 

Bisher ging man bei derartigen Interpplationsschaltungen von der Vorstellung aus, daB fur hochwertige Systeme 
mdglichst ein hochwertiger TiefpaB als analoger Resampler R verwendet werden muB, der alle Signalantelle oberhalb 
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des Nutzsignalbandes B untendruckt. Nach der Erf ihdung werden j edoch fur die Rlterung der digitalen Signale im zweiten 
Interpolationsf liter 4, das in der Folge vereinfacht al's Interpplationsfiiterbezeichnetwird, Filterver^hren hercingezogen, 
die sich ihri wesentlichen kammfilterartig riur auf die Sterbereiche beschrSnken. MaBnahmen der Filterverbesserung, 
die entsprechend der Erfindung diesen Weg beschreiten, haben eine sehr hohe Effektivitat bei der Verbesserung der 

5 Gesamteigenschaften des Systems zur Fotge. Kleinere EinbuGen im Ubertragungsverhalten des Nutzsigraibandes sind 
dabei im Vergleich zuden Verbesserungen im Signal/Rauschverhditnis vernachlassigbar. Auch ist es fur das erfindungs- 
gemaBe Verfahren nicht wichtig, daB die ursprunglichen Abtastwerte auch nach der Interpolation exakt Teil des Signal- 
verlaufs y(t) sind. wie dies beispielsweise bei dem Lagrahge-Verfahren der Fall ist. 

lnRg.4sindschematisch im Zeitbereichdielmpu}santworten h(K=1) bzw. h(K=2) derAbtasthaltefunktionenerster 

10 und zweiter Ordnung K = 1 bzw. K = 2 dargestellt, die der bekannten Abtasthalteschaltung bzw. dem bekannten linearen 
Interpolator entsprechen. Zusatziich sind die Impulsantworten h(K=3) bzw. h(K=4) fur Abtasthalteschaltung dritter bzw. 
vierter Ordnung K=3 bzw. K=4 dargestellt, mit denen das erflndungsgemSiBe Verfahren realisiert werden kann. Aus der 
Darstellung geht auch hervor. daB je, nach der Ordung K der Abtasthaltefunkbbn das VerldnQpfungsfenster ein bzw. zwei 
bzw. drei bzw. vier Abtastwerte des ersfen Abtastrasters erfaBti die im Interpolator 4 gespeichert werden mQssen. Die 

15 zugehorigeh Gewlchtungsfaktoren werden aus dern Zeitdifferenzwert dt und dem Schnittpiinkt mit der jeweiligen Inripuls- 
antwort ermittelt. 

pie jmpulsantwort einer Abtasthalteschalturig der n-ten Ordung hn(t) ist durch folgende Fprmej definiert: 



25 dabei bisdeuten':. 

"n* in hn(t) die Ordnung der rinittigen Nullstellen der zugehorigen Oberfragungsfunktion Hri(0 in jedem Storbereich 
1S.2S. ....pS.... . ' ' ^ 

- "k" ein lalifenderZahlenwert fur die Addition, der die Werte Null bis n annim 
30 "(||)!*stellt deh zugehdrigen Binoniinalkoe^^ 

vet)" bzw. "crCtrkTI) " steht fQr einen Einheitssprung zum Zeitpunkt t bzw. einen um kTakte verschobenen Einheits- 
sprung zurri Zeltpunkl t-kT1 . - . . - : 

Der entisprechende Frequehzyerlauf einiger Ubertrdgurigsfunktionen H(ksn) ist in Fig. 5 dargestellt. Mit dem Para- 
35 meter Kwird der Freqiienzverlauf fur eine Abtasthalteischaltung mit 1 bis 4 Nullsteiien daurgestetlt. J&grOBer die Anzahl 
Kder Nullstellen ist, desto grdBer wird die Dampfung in den Sperrbereichen. Bei der Betrachtung eines vorgegebenen 
Signal/Rauschverhaltnisses. z.B. bei 80 dB, ist erkennbar, daB sich mit steigender Anzahl K an Nullstellen die jeweiligen 
. Sperrbereiche verbreitern. Unabhangig davon nehmen bei festem K-Wert die Sperrbereiche in ihrer Breite mit steigen- 
den Frequen^en zu. Am kritischsten furdas Stdryerhalten sind di(e Storbereiche bei niedereri Frequen^en. Diese kanrien ^ 
. 40 ' jeddch riach dem vort^eschriebeneh Verfahren gezielt mit zusatilich^h Nullstefijen versehen werden. 

Nach dem erf IndungsgeniaBen Verfahren ist-es m6glich, jed§ Impulsaritwdrt h(t) so zu verandern, daB mindestens 
eine weitere Nullstelle durch eine Zeitverschiebung gezielt positioniert werderi kanh. Aus dem Verschiebungssatz d^r 
Fouriertransfbrniation ergibt sich 

-45 f(t -ta)o=> e'^""^^ • F(cd) 

und durch Addition mit der ursprunglichen Funktion f(t) folgt: 

f(t) + f(t - ta) o=> 2e '^"^^^ ■ cos{cDta^) • F( <D ) 

so 

Diese Formel enthalt eine cosinus- Funktion cos( ta • © / 2 ) , dieperiodische Nullstellen aufweist. Da in der Formel der 
fouriertransformierte Signalverlauf F(a)) mit der cosinus-Funktion cbs(ta~- a) /2) multipliziert wird, weist das Produkt 
sowohl die Nullstellen von F (o) als auch diejenigen der cosinus-Funktion auf. Durch die Vorgabe der Zeityerschiebung 
ta kann die erste Nullstelle der cosinus-Funktion so gelegt werden, daB o = cog^^^ ist. Dies ist der Fa|lfur ta = 1 /2f g^^^.. 
55 Die neue Impulsantwort h(t). die von der ursptunglichen Impulsantwort h^itCt) ausgeht, erhSIt man nach dieser 
Betrachtung durch die additive Kbmbination der beiden zeitverschobenen ursprOnglichen Impulsantworten nach fblgen- 
der Beziehung: 

h(t)=h„3^(t) = h^,(t) + h^,(t-ta). 
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Nach diesem Verfahren wird die resultierende Impulsantwort h(t) = h2la(t) von Fig. 6 aus einer Uberlagerung von 
zwei zeitverschobenen Impulsantworten h2i(t) und h22(t) eines tinearen Interpolators h(K:=2) gebildet. Das angewandte 
Uberlagerungsverfahren ist erne Addition der beiden ursprOnglichen Impulsantworten zur neuen Impulsantwort h2ta(t). 
Das Abtastfenster der neuen Furiktion ist grGBer geworden und kann bis zu drei Abtastwerte der ersten Abtastfolge x(k) 
5 . erfassen. SelbstyerstSndlich kQnnen beliebige Impulsantworten h(t) der gleichen Operation untenvbrfen werden, um 
gezielt weitere Nullstellen zu positionieren. So kann beispielswelse die resultierende Impulsantwort von Fig. 6 ihrerseits 
wieder mit einer gleichen. zeitverschobenen Impulsantwort kombiniert werden, um eine zweite zusdtzliche Nullstetle 
gezielt zu positionieren. 

Zum Beispiel sinddrei Obertragungsfunktionen einer Abtasthalteschaltungfurverschiedene Zeitverschiebungen in 
10 Fig. 7 dargestellt. Als Referenz dient die einfache Abtasthalteschaltung erster Ordung Hsh- Nach einer ersten zeitlichen 

Verschiebung um ta = ta1 ergit>t sich der mit HsH+tai dargestellte Kurvenverlaiif. Nach einer zweiten zeitjichen Ver- 

schiebung um ta = ta2 ergibt sich der mit HsH+tai+ta2 dargestellte Kurvenverlauf. 

Beispielsweise wird mit einer ersten zeitlichen Verschiebung von ta = T1 / 2 die Anzahl aller Nullstellen der unge- 

raden StO/stellen veidoppelt Diese StOrstellen enthalten jetzt jeweils zwei Nullstelleii anstatt der ursprOnglichen einen 
IS Nullstelle. Bei einer nochmailgen Anwendung der zeitlichen Verschiebung, aber nun mit ta = 11 /4 wird jeweils eine 

weitere Nullstelle bei bestimmten periodischen StOrstellen, die bei deh geradzahligen Frequenzvielfachwerten 2f 1 , 6f 1 . 

1011. ... liegen, ergdnzt - jeder dieser StOrbereiche weist damit ebenfells zwei Nullstellen anstatt der ursprOnglichen 

einen Nullstelle aul 

Das dritte. oben genannte Vertehren c) zur UnterdrOckung der StOrbereiche ergibt sich, werin in jedem StOrbereich 
20 einzelne Nullstellen jiebeneinander, insbesondere in Squldistanter Anordnung, durch eine Verschiebung irn Frequenz- 
bereich ppsitioniert werden: Dadurch wird die relative Breite der Ddrnpfungszonen.bei den Frequenzvielfachwerten yon 
f1 vergrOBert. wodufch die erforderliche Oberabtastung reduziert werden kann.. Fig. 13 zeigt hierzu im-Firequenzdia- 
gramm ein Beispiel. Die Obertragungsfunktioh Hi((D) einer Abtasthaiteschaltung erster Ordnung wird dabei so modifi- 
ziert. da3 eine heue Ubertragungsfunktion Hneu(^) Produkt^voh zwei im Frequenzbereich verschobenen 
25 ursprOnglichen Obertragungsfunktionen Hait(o>) gebildet wird: 

H„3>) = H2(G>) = HI CO . A© ) • HI ^( CD + Ab> ) 

Dadurch wird das Ddmpfungsverhalten In sdmtlicheri StOrbereichen so ge^lndert. da3 die jeweiligen Ddmpfungszonen 
bei den Frequenzvielfachwerten von f1 relativ breit sind. Obwohl das Durchla(3verhalten yon H2(f) in Fig. 8 bei den 
Frequenzvielfachwerten von f1 ein relatives Maximum aufweist, wird die StOrsignalunterdruckung im ersten StOibereich 
nie schlechter als etwa - 40 dB. Im zweiten StOrbereich liegt die StOrsignalunterdrQckung bereits unter - 50 dB. Soll die 
Stdrsignalunterdruckung weiter verbessert werden,. dann ist es erforderlich die Ubertragungsfunktion H(f) s H3(f) so 
zu modif izieren, daB jeweils eine dritte Nullstelle genau in der jeweiligen StOrbereichsmitte positioniert wird: 

H3(<d) = H1{a)) • H2(a)) . 

Gegebenenfalls werden eine vierte und funfte Nullstelle zwischen die vorhandehen Nullstellen gesetzt und so weiter. 
Dies ist mit dem nachfolgend beschriebeneh Verfahren mdglich. Die resultierende Impulsantwort h2(t) wird durch eine 
Fourier-ROcktransformation wie fotgt gebildet: 

*'• . ■ ■ ■■ ' ' ■' • • ' ■ • ' ' ' ' 

h2(t) = sin(tajl) / col • a(t) + {sin[(2T1-t) • col] - sin(ta)1)} / col - cT(t-tl) + sin[(2T1-t) - ©1] / col • a{t-2T1) 

Die zugehOrige Impulsantwort h2(t) zeigt Fig. 9. Sie hat eine groBe Ahnlichkeit mit der Impulsantwort h(K=2) des linearen 
Interpolators - sie hat etwa die Form eines gleichschenkligen Dreiecks. aber die beiden Schenkel sind. anders als beim 
tinearen Interpolator h(K=2) . leicht konvex gekrOmmt. Dies ist durch die oben beschriebene sinus-Funktion der Impuls- 
antwort bedirigt. 

In Fig, 10 ist scheniatisch als Blockschaltbild eine Anordnung zur Realisierung des Interpolationsfilters 4 dargestellt. 
Die Anordnung entspricht bekannten Transversalfifteranordnungen mit zeitabhangigen Kbeffizienten. Sie enthait eine 
Speichereinrichtung rinr. die eine bestimmte Anzahl 1 bis L bzw. kn bis kn-L+1 Abtastwerte des ersten Abtastrasters x(k) 
speichert und diese Werte mittels einer Gewichtungseinrichtung g mit Gewichtungsfaktoren Co bis Cl.^ multipliziert und 
in einer Addiereinrichtung ad aufaddiert. 

Die Bildung der jeweiligen Gewichtungsfaktoren ergibt sich aus der zeitverschobenen Irhpulsantwbrt nach fblgender 
Gleichung: 

C ,(dt) = h(l-T1 + dt) mit I = 0. 1 , 2. .... L-1 

Im Unterschied zu den bekannten Verfahren werden nach der Erfindung die zeitabhangigen Gewichtungsfaktoren (= 
Koeffizienten) nicht einer gespeicherten Tabelle entnommen oder mittels einer umfangreichen Berechnung ermittelt. 
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sondern auf einfache Weise mittels einer Matrix mx gebildet deren Schaltung in Fig. 1 3 ausfQhrlich dargestellt ist. Diese 
Realisierung ist nur fur die Impulsantworten nach den erfindungsgemdlBen Veiiahren a) und b) mdglich, weil diese 
abschnittsweise.polynominelle Funktionen llefern. 

Der Zeitdifferenzwert dt wird einer Matrixstufe mx zugefuhrt, in der alle Komponenten def Gewichtungsfaktoren Ci 
5 bzw. C| matrixfarmig gebildet werden. Die einzelnen Matrixwerte, z'B die Konstanten bq q bis a^,^ i_ -i Fig. 13. werden 
aus einer Speichereinrichtung ausgelesen und mit dem jeweiligen Zeitdifferenzwert dt bzw. mit der entsprecherxien 
Potenz von dt, dt^, dt^. multipliziert und aufaddiert und bilden schlie6lich die zatabh§ngigen Gewichtungsfaktoren 
Co(dt)bisCL.i(dt). 

In Fig. 11 ist eine besonders vorteilhafte Schaltungsanordhung fur das Interpolatpnsfilter 4 dargestellt, bei dem 
70 zusatzlich eine Umschaiteinrichtung ux und eine erweiterte Matrixstufe mxx vorhanden ist. Dieses Ausfuhrungsbeispiel 
des Interpolationsf liters 4 erweist sich dann von Vorteil, wenn die Ubertragungsfunktion H(f) eine zusatzliche Nullstelle 
' nach dem zweiten Verfahren b erzielen soil. Zum Belspiel lessen sich die beiden Hnearen Abschnitte des linearen Inter- 
polators abschnrttsweise jeweils als Gerade beschreben: Dies erfelchtert den Aufwand fur die Berechnung der jeweiligen 
Gewichtungsfaktoren. In Fig. 11 enthdlt die Matrixstufe mx beispielsweise die beiden Matrizen A' und A", dehen die, 
IS Abscbnittsbereiiche dt < ta und dt s ta entsprechen. In Fig. 6 ist beispielsweise eine Impulsantwort h2ta(t) dargestellt, 
die abschnittsweise aus geraden StOcken zusammehgesetzt ist. Die gespeicherten Abtastwerte in Fig. i1 umfassen 
L+1 Abtastwerte von 1 bis L+1 . Ihnen entsprechen die gespeicherten Abtastwerte x(kn) bis x(kn-L). . ^ 

Am Beispiel der modiflzierten Impuisanwort h2ta(t) von Fig. 6 soil die Bestimmung der Gewichtungsfaktoren Co. 
Ci. O2 entsprechend.dem Zeitdifferenzwert dt erldutert werden, wobej die Berechnung. abschnittsweise erfolgen soil. 
20 Dazu wird die Impulsantwort h2ta(t) in'fQnf Abschnitte zerlegt. AuBerhalb dieser Atechnitte hat die Impulsantwort den 
WertNuil; V \ . ^ : - ; 3 V 

- - AbsdinmyonO<t<ta: h(t) = t ^ ^ i . . i. 

2S . Abschnittvonta <t < 1 : h(t) = -ta + 2t 2. 

Abschnittvon 1 <t < 1+ta : , h(t)=2-ta 3. 

Abschhittvon1+ta<t<2: h(t)=4 + ta-2t 4. 

Abschnittvon2<t<2+ta: h(t)-2 + ta-t .5. 

For die nachfolgende Berechnung der Gewichtungsfaktoren gelten tolgende Gleichungen: - 

-35 C,(dt) = h(l+dt) mit dem laufenderi Parameter 1 = 0. 1.2, ...und 

h(t) = h^^(t) + h^(t-ta) : 

. Dielrr43ulsaritworth2ta(t) hat beiden Frk|uenzwerten pf1. 1f1 und 2f1 jeweils eine Uhstetigkeitsstejle. Da auch bei ^ 
'40 ZeitMrerten 0+ta. T1+ta urid 2T1fta jeweils eine Unstetigkeit vorliegt ist die rnathematische Besdhreibung der geraden 
Zwischenstucke bjs zurh jeweiligen Zeitwert ta abschnittsweise viel eirifacher,' Dies zeigen die bbigen funf einfachen^ 
Impulsantworten h(t). Sie liefeirn folgehde Gewichtungsfaktoren fur die funf mdglichen Zeltbeireiche von dt innerhalb der 
Inrpulsantwort: ^ ' . 



45 



so 



0 < dt < ta bzw. ta <dt< 1 

C,(dt)= dt C,(dt) = -ta + 2dt 

C,{dt) = 2.ta C,(dt)= 2 + ta-2dt 

C,(dt)= ta-dt C,(dt)= 0 



55 
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Dies ergibt fOr die Matizen A' und A" folgende umschaltbare Koeffizienten oder Matrixwerte, die als Faktoren fOr die 
nachfolgende Multlpllkation bzw. als zeitunabhdngiger Wert einer Speichereinrichtung entnommen werden: 



5 




0 


1 




-ta 


2 






2-ta 


0 


A" = 


2+ta 


-2 






ta 


-1 




0 


0 
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Fur das dritte Verfahren c) mit der Verschlebung im Frequenzbereich ist die Berechnung der Gewichtungsfakoren nicht 
mehr mittels einer Matrixanordung wie in Fig. 10 Oder Fig. 1 1 mdgiich. sondern sie erfordert eine Berechung von sinus* 
Werten oder eine Entnahme dieser Werte aiis einer gespeicherten sinus-Tabelle. Denn die erfbrderlichen sinus-Werte 

15 lassen sich nicht durch eine einfache Polynomreihe'approximieren - in diesem Falle wQrde man direkt die Polynomreihe 
und damit die MatrixlOsung verwenden. 

. In Fig. 1 2 ist etwas ausfuhrlicher als in Fig. 2 dargestellt, wie aus einer ersten Abtastfolge x(k) eine zweite Abtastfolge 
y(d) mit beliebigem Frequenzvertidltnis ierzeugt wird. Der gesuchte Abtastwert Pd, der zu dem Interpolations- Oder Abt- 
astzejtpunl^ tP gehOrt wind.durch eine. insbesondere polynbmine|le, Interpolation aus der erstep Abtastfolge x(k) gebil- 

20 det. Hierzu wird der Zertdifferenzwert dt aus dem nachstliegenden Abtastwert Pn der ersten Abtastfolge x(k) und dem 
gewdnschten Abtastzeitpunkt tP bestimmt. Entsprechend diesem Zeitdifferehzwert tP vyerden die Gewichtungsfaktpren 
Ci aus'der Obertragungsfunktion h(t) ermittelt. Durch die Speichereinrichtung mr wird dabei sichergestellt, daB die ^Ipr- ■ 
derlichen Abtastwerte Ph und Pn+1 fOr die Interpolatiqnsberechnung von Pd zur VerfOgung stehen^ Diei Berechnurig - 
der vorausgehehden und nachfdgenden Abtastwerte Pd-2 bis Pd-i-2 erfolgt auf dhnfiche Weise durch VerschietMjng des 

25 Interpolationsfensters und Bestimmung des jeweillgen Zeitdiffererzwertes dt. 

PatehtansprQche 

1. Verfahren zur digitalen Interpolation von Signalen durch Multiplikation verzdgerter. mit einem ersten Taktsignal (f 1) 
30 verkoppeiter Eingangswerte (x(k)) mit ehtsprechenden Gewichtungsfaktoren (Ci; C)) und Aufaddierung der muttipli- 
zierten vyerte, die von einem Zeitdifferenzwert (At; dt) abhangig sind, der durch den tnterpolatibnszeitpunkt (tP) und 
das Zeitraster des ersten Taktsignals (f1) bestimmt wird, 
aekennzelchnet durch fblo ende Merkmale: 

35 - mittels einer Impulsantwort h{t) im Zeitbereich werden die Gewichtungsfaktbren (Ci; Cj) bestimmt, wobei die 
Obertragungsfunktion H(f ) im Frequenzbereich ein Signaldampfungsverhalten aufweist. das bezuglich der 
Sperrbereiche im wesentlichen auf die bei den Frequenzvielfachen des. ersten Taktsignals (f 1) liegenden StOr- . 
bereiche (1S, 2S, ... pS. ...) beschrankl ist. ! 
jedem dieser Storbereiche werden mindestens zwei Nullstellen der Obertragungsfunktion H(f) im Frequenzbe- 

40 reich zugeordnet, wobei r 

die mindestens zwei Nullstellen in jedem Stdrbereich nebeneinander liegen o^ 

mindestens einern der Stdrbereiche (pS) mindestens eine weitere Nullstelie der Obertragungsfunktion H(f) 
zugebrdnet wird, yyenn die zwei Nullstellen als Nullstelie zweiter Ordnung (K=2) aufeinand^ liegen. 

45 2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, daB eine neue Impulsantwort h(t) , deren erste weitere Nullst- 
elle der Obertragungsfunktion H(f) im Frequenzbereich bei einer vorgegebenen Frequenz fn liegen soli, nach fol- 
gender Bemessungsregel aus einer ursprunglichen Impulsaritwort hait(t) 

- h(t) =-h W = h attW + h aii{t - ta).mit ta = 1 / 2fn gebildet wird. 
so • ' . 

3. Verfahren nach Anspruch 2, dadurch gekennzeichnet. daB zusStzliche weitere Nullstellen der Obertragungsfunktion 
H(f) durch eine wiederholte Anwendung des Verfahrens gemdB Anspruch 2 beliebig positioniert werden kOnnen. 

4. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB jedem Stdrbereich (1S, 2S pS, ...) mindestens eine 

55 weitere Nullstelie der Obertragungsfunktion H(1) im Frequenzbereich zugeordnet ist, so daB alle Stdrberelche min- 
destens drei Nullstellen aufweisen. 

5. Verfahren nach Anspruch 4, dadurch gekennzeichnet. daf3 die mindestens drei Nullstellen mittig zu den jeweiligen 
StCrbereichen (1S, pS, ...) bei den Frequenzvielfachwerten (1f1 , 2f1 ..... pf1, ...) des ersten Taktsignals (f1) liegen. 
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6. Verfahren nach Anspruch 5, dadurch gekennzeichnet. daB die Impulsantwort h(t) einer Abtaslhaftefunktion hn(t). 
bzw. h(K=n) mit mittig liegeriden Nullstellen. die mindestehs die Ordnung K = n = 3 aufweisen, durch folgende. 
Forme! definiert yvird: . . 



^^^^ ' {n-iyi ^ ^ (*) (t-Jtri)«-^ X a{t-kTl)] 



dabei bedeuten: 

"n" die Ordung der mittigen Nullstellen der Ubertragungsfunktion Hn(f). - . 

"k" ein laufender Zahlenwert fur die Addition-, der die Werte Null bis n annimmt urid 

V(t)" bzw. "cT(t - kTI) " einen.Einheitssprung zum Zertpunkt t bzw. zum verzdgerten Zeitpunkt t - kTI. 

7. Verfahren nach Anspruch 6, dadurch gekennzeichnet,. da B aus einer alten Impulsantwort hnajt(t) eihe neue Impids- 
antwort hnneu(Q 961^^3 einem der Verlahren von Anspruch 2 Oder 3 gebiM^^ 

8. Verfahren nach Anspruch 1 , dadurch gekennzeichnet; daB aus einer ursprunglichen Ubertragungsfunktion Mah(f) 
durch eine Verschiebung im Frequenzbereich um die Frequenzwerte + Af und - Af eine.neue Clbertragungsfuhktion 
HneuCP' n^ch^folgender Formel . ; ; v . . P ' 

gekDildet wird, wpdurch derjeweilige Sperrbereich gegenOber einer rnittigen Nullstellenlage verbreitert wird: 

9. Verfahren nach einem der Anspruche 1 bis 7, dadurch gjBkennzeichnet da6 zusStzliche Nullstellen der Ubertra- 
gungsfunktion. H(f) durch Wiederholung des Verkhrehsger^ 

10. Verfahren hach einem der Ansp'rOche 1 bis 7, dadurch gekennzeichnet daBdie Gewichtungs^ktoren (Co ..... C|_. 
i)mittels eirierMatrixstufe (nrix: rnxx) gebildet werden. 

11. Verfahren nach einem der AnsprQche 1 bis'7 oder Anspruch 10, dadurch gekenrizeichnet, daB die Impuisaritwort 
h(t) fOr die Berechnunig der Gewichtungsfaktoren (Ci; C|; Co. C^.i) abschnittsweise bezQgtich dem Zeitdifferenz- 
wert (dt) def iniert ist. ' 

12. Verfahren nach eihem der Anspruche 1 bis 7 und/oder einem der Anspruche lOoder 11 , dadurch gekennzeichnet, 
daB die imjsulsantwort h(t) rriit einer FIR-Filterstruktur realisiertvWrd. wobei die GewichtungsfiaKtoren mittels etner 
lyiatrixstufe (mx; n^} ge^^ • . 

13. Verfahren nach Anspruch 12, dadurch gekennzeichnet. daB die Gewichtungsfaktdr bei einer abschriittsweise 
def inierten impulsantwort .h(t) mittels einer umschd^ 
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